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Description 

Parametric recursive digital filter 



The invention relates to a parametric recursive digital 
filter having at least one delay unit and a coefficient 
network connected to the delay unit or the delay units, 
respectively. 
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10 Digital filters having different filter structures and 
in various configurations are known, for example, from 
U. Tietze, Ch. Schenk, Halbleiterschaltungstechnik, 
(semiconductor circuit technology) , 9th edition, 1991, 
pages 791 to 851, and consist of at least one delay 

15 unit formed of a delay element and one coefficient net- 
work determining the filter structure and the transfer 
characteristic including cut-off and centre frequen- 
cies. In numerous applications, it is necessary to 
change the cut-off or centre frequency of a filter in 

20 dependence on a control quantity (parametric filter) . 

In conventional digital filters, the cut-off or centre 
frequency is controlled either with the aid of coeffi- 
cient tables or real-time calculation of the filter co- 
25 efficients. This results in very high computing or mem- 
ory complexity, particularly in the higher-order fil- 
ters . 



In W. Schussler, W. Winkelnkemper , "Variable Digital 
30 Filters", Arch. Elektr. Ubertr. , Vol. 24, 1970, is- 
sue 11, pages 524-525, an arrangement is proposed in 
which the delay elements of a conventional digital fil- 
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ter are replaced by all -pass filters in order to obtain 
a parametric filter. 

At the end of the article, however, the authors concede 
5 that, as a rule, this procedure cannot be implemented 
in digital filters since the all-pass filters comprise 
an attenuation when the highest all-pass filter coeffi- 
cient is not equal to zero. As a solution for this 
problem, it is suggested either to use only all -pass 

10 filters in which the highest all-pass filter coeffi- 
cient is equal to zero or only to use non-recursive 
filter structures. Using this proposed solution, how- 
ever, only a few and only special filters can be imple- 
mented so that the possible applications are very lim- 

15 ited. 

It is the object of the invention to specify a paramet- 
ric recursive digital filter which can be implemented 
without the abovementioned restrictions. 

20 

The object is achieved by a parametric recursive digi- 
tal filter according to Patent Claim 1. Embodiments and 
further developments of the concept of the invention 
are the subject-matter of subclaims. 

25 

It is an advantage of the invention that the cut-off 
frequencies of recursive digital filters of otherwise 
any other configuration can be controlled by changing a 
single coefficient and this one coefficient can be de- 
3 0 termined by a very simple method of calculation. In 
particular, a very uncomplicated calculation is ob- 
tained from the fact that, as an approximation, a lin- 
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ear relationship exists between this one coefficient 
and the cut-off or centre frequency, respectively. 

The parametric digital filter according to the inven- 
5 tion can be implemented by means of any recursive fil- 
ter structure including digital wave filter structures. 
It is also advantageous that in each case a filter ele- 
ment with controllable phase angle (e.g. all-pass fil- 
ter, Laguerre-Kautz-Filter) is in each provided as the 

10 delay unit. The filter structure of the digital filter 
according to the invention is determined by the coeffi- 
cient network as is customary, but the associated cut- 
off f requency/centre frequency can be adjusted via the 
phase angle of the filter element or of the filter ele- 

15 ments, that is to say preferably via only a single co- 
efficient . 

According to the invention, in detail, a parametric re- 
cursive digital filter of the type initially mentioned 
2 0 is proposed which comprises 

a delay unit, the phase angle of which is control- 
lable and which comprises a delay element and a 
phase network interconnected therewith, 
25 - a positive feedback network connected to the delay 
unit and compriseing a positive feedback path, 
a feedback network connected to the delay unit and 
compriseing a feedback path, 

the feedback path of the feedback network being 
30 connected to the output of the delay element in the 

delay unit, and 
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the cut-off frequency/ centre frequency of the fil- 
ter being adjustable via the phase angle of the de- 
lay unit . 

5 In more complex filters and in higher-order filters, in 
particular, a plurality of delay units is used and the 
delay units can be formed preferably by identical fil- 
ter elements selected in the same manner. As a result, 
all filter elements have the same coefficient so that 
10 the calculation of the coefficients only needs to be 

performed once and can be used for all filter elements. 
This saves considerable computing expenditure. 

The delay units can comprise a plurality of, in par- 
15 ticular two, delay elements. In addition, the positive 
feedback network can comprise a plurality of positive 
feedback paths and the feedback network can comprise a 
plurality of feedback paths. 

20 In each case, an all-pass filter is preferably provided 
as a delay unit, and frequency-influencing filter units 
such as low-pass, bandpass, high-pass filters etc. can 
also be used with corresponding restrictions and/or 
special features in the characteristics. 

25 

An all-pass filter preferred for use as delay unit com- 
prises the following: 

a first adder, one input of which forms the input 
30 of the delay unit, 

a second adder, the output of which forms the out- 
put of the delay unit, 
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a coefficient section which is connected between 
the output of the first adder and a first input of 
the second adder, 

a first delay element which is connected between 
the input of the delay unit and a second input of 
the second adder, 

a second delay element which is connected between 
the output of the delay unit and a second input of 
the first adder, 

the phase angle of the filter element being adjust- 
able by changing the coefficient of the coefficient 
section, and 

the output of the first and/or second delay element 
being provided for connecting a feedback path. 

Apart from this embodiment of an all -pass filter which 
is of advantage for many applications, however, any 
other embodiment of an all-pass filter can also be 
used . 

To achieve a low- expenditure filter, two delay units 
comprising delay elements can be interconnected with 
one another (in each case) in such a manner that only a 
total of three delay elements are provided, one delay 
element being used for both delay units. 

In the text which follows, the invention will be ex- 
plained in greater detail with reference to the illus- 
trative embodiments shown in the figures of the draw- 
3 0 ing, in which: 

Figure 1 shows the signal flowcharts (structure) of 
two first-order all-pass filters, 
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Figure 2 shows the phase response or frequency bend- 
ing function, respectively, of a first-order 
all-pass filter as a function of the all- 
pass filter coefficient, 



10 



Figure 3 shows the frequency resolution of the para- 
metric filter for various all-pass filter 
parameters , 

Figure 4 shows the structure of a conventional finite 
impulse response (FIR) filter, 



15 



Figure 5 shows the structure of a corresponding para- 
metric finite impulse response filter, 



Figure 6 shows the detailed, computationally effi- 
cient structure of a parametric FIR filter, 

20 Figure 7 shows the structure of a biquad filter as an 

example of a second-order infinite impulse 
response (IIR) filter before mapping, 

Figure 8 shows the structure of a biquad filter as an 
25 example of a second-order infinite impulse 

response (IIR) filter after mapping, 



30 



Figure 9 shows the structure of parametric infinite 
impulse response filters with a low-pass 
filter as the first delay unit, 



Figure 10 shows the pole -zero shift in the Z domain as 
a function of the all-pass filter parameter 
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in a second-order parametric IIR filter 
structure according to the invention, 



Figure 11 shows the pole- zero diagram of a prototype 
5 filter according to the invention, 

Figure 12 shows the amplitude /frequency response of 

the prototype filter according to Figure 11 , 
and 

10 

Figure 13 shows a computationally efficient version of 
a second-order parametric IIR filter accord- 
ing to the invention, 

15 Figure 14 shows the impulse response of a target fil- 
ter for a second-order parametric IIR filter 
according to the invention, 

Figure 15 shows the resultant amplitude/ frequency re- 
20 sponse of the target filter, 

Figure 16 shows the pole- zero diagram of the target 
filter. 



25 Digital filters in which simple delay sections are re- 
placed by frequency- dependent delay sections such as, 
e.g. a first-order all-pass filter, are often called 
frequency distorting filters since they represent fil- 
ters, the specifications of which are implemented on a 

30 distorted nonuniform frequency axis. 
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A shift in the amplitude/ frequency response |H(z) | can 
be achieved if z is replaced or mapped by an expression 
W(z), where W(z) must meet the following conditions: 

5 - the inside of the unit circle (Z plane) must be 
mapped again to the inside of the inner circle, 
the unit circle must be mapped to itself. 

The first condition means that, if H(z) is stable, 
10 H(W(z)) is also stable. The second condition makes it 
possible to shift or to "map" the amplitude/frequency 
response into another arbitrary frequency range. 

W(z) must be clearly an all-pass filter since this 
15 meets both the above conditions. The amplitude/ fre- 
quency response |H(z) | of the all-pass filter which is 
constantly equal to one as is known ensures that the 
amplitude/frequency response of the overall filter is 
not changed and that the first condition is met, the 
20 phase response effecting the actual frequency displace- 
ment, i.e. being responsible for meeting the second 
condition. Thus, a displacement or bending of the fre- 
quency axis occurs, the displacement of which can be 
controlled by the filter coefficient of the first-order 
25 all-pass filter which is why the filter efficient y of 
the first-order all-pass filter is also called fre- 
quency curvature parameter. 

Figure 1 shows two examples of first-order all -pass 
30 filters. 

Example (a) comprises a first adder 1, one input of 
which forms the input of a delay unit 2, and a second 
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adder 3, the output of which forms the output of the 
delay unit 2. A coefficient section 4 is connected be- 
tween the output of the first adder 1 and a first input 
of the second adder 3. A first delay element 5 is con- 
5 nected between the input of the delay unit 2 and a sec- 
ond input of the second adder 3. Finally, a second de- 
lay element 6 is also connected between the output of 
the delay unit 2 and a second input of the first ad- 
der 1 . 

10 

The second example (b) again comprises a first adder 7, 
one input of which forms the input of a delay unit 9, 
and a second adder 8, the output of which forms the 
output of the delay unit 9* A first coefficient sec- 

15 tion 10 is connected between the output of the first 

adder 7 and a first input of the second adder 8. A sin- 
gle delay element 11 is connected between the output of 
the first adder 7 and a second input of the second ad- 
der 8. In addition, the output of the delay element 11 

20 is conducted to the second input of the first adder 7 
via a second coefficient section 12. The first coeffi- 
cient section 10 comprises the coefficient -y and the 
second coefficient section 12 comprises the coefficient 

y • 

25 

A first-order all-pass filter has the following trans- 
fer function: 

\-y*z' 1 
z -y 

30 The first-order all-pass filter changes its phase re- 
sponse as a function of the all-pass parameter y which 
can move within a range of -1 < y < +1. The change of 
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this frequency response, in turn, corresponds to the 
bending of the frequency axis or, respectively, its 
mapping into a new frequency range. 

5 Figure 2 shows the frequency response or frequency 

bending function of the first-order all -pass filter as 
a function of the all-pass coefficient or frequency 
curvature parameter y. It can be seen that high fre- 
quencies are mapped to low frequencies for negative y 
10 values and low frequencies are mapped to high frequen- 
cies for positive y values. The greater the amount of y, 
the greater the frequency bending. There is no fre- 
quency bending at y = 0, i.e. the frequency axis is 
mapped to itself again in this case (bisector) . 

15 

Since delay elements in traditional filter structures 
do not appear in the form of z but in the form of z" 1 , 
it is also appropriate for filters with frequency bend- 
ing to work not with W(z) but with W _1 (z) . In this way, 
20 it is impossible to replace all delay elements z" 1 with 
their frequency- dependent counterpart 

W'(z) = D(z)= 

25 in any filter structures. Although it is not quite cor- 
rect that this applies to any filter structures, this 
question and its significance will be discussed in de- 
tail later. 

3 0 Due to the bending of the frequency axis, some fre- 
quency ranges are resolved more finely and others more 
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coarsely in dependence on the frequency curvature pa- 
rameter y. 



If a filter has a local frequency resolution of Af (f ) 
5 before the "bending", the bent filter has a local fre- 
quency resolution of Af w (f ) / where the following holds 
true : 



10 



a/ w (/) = Af(f w (f,r))* 1 7 \+ n + f ifrt 

fs 



where f s is the sampling frequency and f w (f / y) stands 
for a function which specifies the corresponding fre- 
quency in the bent frequency range for. a given fre- 
quency and all -pass parameter. The frequency of the pa- 
is rametric filter as a function of the frequency and of 
the all-pass parameter y is: 



y*sm( — ) 

/„(/» = / + * arctan( /' ) 

JS 



2 0 The turning point frequency f tp is the frequency at 

which the frequency resolution of the parametric filter 
is the same as before the frequency bending. It can be 
shown that the following holds true for f tp : 



25 



f = * arccos(y) 
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For frequencies below f tp , the frequency resolution is 
improved with a positive y. It worsens correspondingly 
for negative y. Figure 1 shows in this respect the fre- 
quency resolution of the parametric filter for various 
5 all -pass parameters y. The turning point frequency f tp 
is given when the frequency resolution is equal to one. 

In an example of a traditional FIR (Finite Impulse Re- 
sponse) filter structure according to Figure 4, two de- 

10 lay elements 13, 14 are connected in series and driven 
with an input signal x [n] , the taps of the delay ele- 
ments 13, 14 etc. being supplied to an adder 18 at the 
output end via coefficient sections 15, 16, 17 etc., 
forming an output signal y [n] . The structure of a fre- 

15 quency-bending FIR filter is obtained by replacing each 
delay element 13, 14 etc. occurring in the traditional 
FIR filter structure (e.g. Figure 4) with a frequency- 
dependent delay element 19, 20 having a transfer func- 
tion D (z) . 

20 

The structure of such a frequency-bending FIR filter is 
shown in Figure 5 . 

If then the frequency-dependent delay sections are re- 
25 placed by their corresponding all-pass filter structure 
(e.g. by the structure shown in Figure la) and the 
structure is simplified by dissolving redundant 
branches, a computationally efficient structure of a 
frequency-bending FIR filter is obtained. 

30 

To this end, Figure 6 shows the detailed computation- 
ally efficient frequency-bending FIR structure for the 
general structure according to Figure 5. This structure 
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comprises series-connected delay elements 21, 22, 23 
etc., an addition node in the form of an adder 24 , 25 
etc- being inserted between in each case two delay ele- 
ments 21, 22, 23 etc. The signals for the coefficient 
5 sections 21, 22, 23 etc. are picked up at the input of 
the delay element 21 and the outputs of the adders 24, 
25 etc. The other input of each of the adders 24, 25 
etc., which is still free in each case, is connected to 
the output of a subtractor 28, 29 etc. via a coeffi- 
10 cient section 26, 27 etc. to which the coefficient y 
(all-pass filter parameter y) is applied. 

The inputs of the subtractor 28 are connected to the 
input of the delay element 21 and to the output of the 
15 delay element 22. The inputs of the subtractor 29 are 
connected to the input of the delay element 22 and to 
the output of the delay element 23. The result is that, 
for example, the delay element 22 is used for two delay 
units (e.g. delay units 19 and 20 from Figure 5) . 

20 

Replacing a delay element in one of the traditional IIR 
filter structures with a frequency-dependent (disper- 
sive) delay element as has already been explained for 
the FIR filter provides a structure which contains 
25 feedback branches which are free of delay but which 
cannot be implemented in this way, as has been ex- 
plained in the prior art mentioned initially. 

Figures 7 and 8 show a recursive biquad filter before 
3 0 and after the mapping. 

Figure 7 is based on the arrangement according to Fig- 
ure 4 due to the fact that the structure has been lim- 
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ited to a second-order structure and has been expanded 
_by a feedback network. The feedback network comprises 
three coefficient sections 30 , 31 and 32 which, on the 
one hand, are connected to the taps of the delay ele- 
5 ments 13, 14 and, on the other hand, to two inputs of 
an adder 33 (coefficient sections 31 and 32) and to the 
output of the adder 33 (coefficient section 30) . A 
third input of the adder 33 is driven with the input 
signal x [n] . 

10 

Correspondingly, the filter shown. in Figure 8 is based 
on the filter according to Figure 5, that is to say by 
limiting the structure to a second-order structure and 
expanding it by the same feedback network as in Figure 
15 7. 

A WIIR filter structure which is also ^suitable for 
higher filter orders provides for all other delay ele- 
ments apart from the first one to be replaced by first - 

20 order all-pass filters,, the first one being replaced by 
a first-order low-pass filter. This measure makes it 
necessary to recalculate (to map) the filter coeffi- 
cients of the original IIR filter to the new structure 
which is done by means of a relatively simple recursive 

25 formula. Figure 9 shows the corresponding new IIR fil- 
ter structure (Laguerre filter) . 

The achievable structure according to the invention now 
consists in replacing all delay elements with all-pass 
30 filters (particularly of the first order) . However, to 
avoid the zero-delay feedbacks, feedback is not ef- 
fected after the respective dispersive delay units (in 
the prior art) but, according to the invention, di- 
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rectly after the delay element within a dispersive de- 
lay unit. The forward branch (FIR part of the filter) 
with its coefficients b 0/ b lt b 2 etc. remains unaffected 
by this measure. 

5 

In the text which follows, the procedure in designing a 
second-order frequency-bending IIR filter is explained 
with the design of a 10 Hz high-pass filter. 

10 Using frequency-bending filters, a fixed ''prototype 
filter" can be designed which has an arbitrary fixed 
cut-off frequency fg and which can then be displaced 
into any desired frequency range with the aid of the 
frequency curvature parameter y. 

15 

Figure 10 shows the pole-zero shift in the Z domain as 
a function of the frequency curvature parameter y. 

Since originally filters are assumed which are designed 
20 in the direct form, it is appropriate to utilize the 

advantages of the direct form as far as possible. As is 
already known, the frequency range within which the di- 
rect form has particularly good characteristics in the 
sense of quantization sensitivity and dynamic range, is 
25 around f a /4 . For this reason, the cut-off frequency of 
the prototype filter is firmly left at f g = f a /4. 

Figures 11 and 12 show the associated pole- zero diagram 
and the amplitude/ frequency response of the prototype 
30 filter. 



Next, the frequency curvature parameter y must be cal- 
culated so that the prototype filter is shifted towards 
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the required cut-off frequency. For this purpose, the 
abovementioned equation 



2*7T* f 

y * sin( — ) 

fs / fs 



f B re(f>Y) = f + — *arctan( ^ fc „ 

l-r*cos( — -^-) 

^5 



can be transformed in such a manner that the necessary 
frequency curvature parameter y can be calculated from 
it. The calculation formula which produces the fre- 
quency curvature parameter y with the aid of the fixed 
cut-off frequency of the prototype filter f pre and the 
desired cut-off frequency f g is: 



f 

Y = -( Is. . ) 

r \ 2*7r* f 2*n*f (f -f)*x' 

sin( + cos( J * ) * tan r 7 '" Jg) ) 

./a ./a J a ' 



If f pre is firmly left at f a /4, the previous equation is 
reduced to: 



4 fa 



The expression f pre for the fixed cut-off frequency of 
the prototype filter was selected since it can be said 
that this filter is designed in the distorted frequency 
range. One can return from the distorted frequency 
range into the original frequency range by operating 
the parametric filter with the frequency curvature pa- 
rameters y calculated above. 



I 

I 
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Figure 13 shows a computationally efficient embodiment 
of a second-order parametric IIR filter. This embodi- 
ment is obtained by expanding the arrangement according 
to Figure 6 correspondingly by one feedback network. 

5 

From the second-order parametric IIR filter, the im- 
pulse response shown in Figure 14 and the resultant am- 
plitude/frequency response of the target filter shown 
in Figure 15 is obtained. The associated pole -zero dia- 
10 gram of the target filter is shown in Figure 16. 

The abovementioned digital filters are preferably im- 
plemented in microprocessors, signal processors, micro- 
controllers, computing devices etc. The individual fil- 
15 ter components such as, for example, delay units, delay 
elements, coefficient sections are then hardware compo- 
nents of the microprocessors, signal processors, micro- 
controllers, computing devices, etc. which are corre- 
spondingly used by means of software. 



20 
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Patent claims 

1. Parametric recursive digital filter comprising a 
delay unit, the phase angle of which is controllable 

5 and which comprises a delay element and a phase network 
interconnected therewith, 

a positive feedback network connected to the delay 
unit and compriseing a positive feedback path, 

a feedback network connected to the delay unit and 
10 compriseing a feedback path, 

the feedback path of the feedback network being 
connected to the output of the delay element in the de- 
lay unit, and 

the cut-off frequency/centre frequency of the fil- 
15 ter being adjustable via the phase angle of the delay 
unit . 

2. Filter according to Claim 1, in which a plurality 
of delay units are provided. 

20 

3. Filter according to Claim 2, in which the delay 
units are identically designed and are controlled in 
the same manner. 

25 4. Filter according to one of the preceding claims, in 
which the delay unit(s) comprise (s) a plurality of de- 
lay elements. 

5. Filter according to one of the preceding claims, in 
3 0 which the positive feedback network comprises a plural- 
ity of positive feedback paths. 
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6. Filter according to one of the preceding claims, in 
which the feedback network comprises a plurality of 
feedback paths. 

5 7. Filter according to one of the preceding claims, in 
which in each case an all-pass filter is provided as 
delay unit (s) . 

8. Filter according to Claim 7, in which at least one 
10 delay unit constructed as all -pass filter comprises the 

following : 

a first adder, one input of which forms the input 
of the delay unit, 

a second adder, the output of which forms the out- 
15 put of the delay unit, 

a coefficient section which is connected between 
the output of the first adder and a first input of the 
second adder, 

a first delay element which is connected between 
20 the input of the delay unit and a second input of the 
second adder, 

a second delay element which is connected between 
the output of the delay unit and a second input of the 
first adder, 

25 the phase angle of the filter element being adjust- 

able by changing the coefficient of the coefficient 
section, and 

the output of the first and/or second delay element 
being provided for connecting a feedback path. 

30 

9. Filter according to one of Claims 2 to 9, in which 
two delay units comprising delay elements are (in each 
case) interconnected with one another in such a manner 
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that only a total of three delay elements are provided, 
one delay element being attributable to both delay 
units . 

10- Filter according to one of Claims 1 to 6, in which 
a frequency- influencing filter unit is provided as de- 
lay unit • 



} 
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Abstract 

Parametric recursive digital filter 

Filter comprising a delay unit, the phase angle of 
which can be controlled and which comprises a delay- 
element and a phase network interconnected therewith, 
and with positive feedback and feedback networks con- 
nected to the delay unit, the feedback path of the 
feedback network being connected to the output of the 
delay element in the delay unit and the cut-off fre- 
quency/centre frequency of the filter being adjustable 
via the phase angle of the delay unit. 

15 'Figure 13 
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a) 



b) 



FIG 1 




Frequency (normalized to Nyquist frequency [fa/2]) 
FIG 2 
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FIG 4 



FIG 5 
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5/7 



P/N Diagram of the prototype filter 




FIG 11 



Amplitude/frequency response of the prototype filter 
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FIG 12 
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Comparison of the ideal impulse response with the 
x 10 -4 real impulse response of the target filter 




FIG 14 



1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 
Samples 





0 




-5 




-10 




-15 


m 


-20 


.E 




o 

"O 


-25 


Z3 

S 




1 


•30 




-35 




-40 




-45 




-50 



Amplitude/frequency response of the target filter 
(calculated from impulse response) 







t i • 1 1 ty i t t t t tin 




i i it it 


r i i 


1 1 it i 






* > I >jy t i t MMMI 














L J iU...l.J.Vlt k|J k . - . 
i/miim I • t • mil 

t/ • » • • 1 » HI 

• nfi-rtirt- * - 1 "*»"•" J [f|]f 

1 1 1 ttM • * 1 « 1 t •<• 
*•••!•• 1 1 1 « 1 1 1*1 

. J -•••J IU..«l.J.llll»Jt 

■ -l-r-rtin- • - 4 - "i -r t r rnr 
» t • • »•»••■• 

1 1 • » III • • I 

1 1 II III * 1 1 1 tllll 


I > I 
k.kl 

• t • 

I- - r i 

• • i 
t i • 


• i »i *• 

• 1 M II 

-X-!r!i — 

^JJIU ... 

_»_»JLU 

-!-!-» rn 

I Mill 


- t - r r 

• 4 .W 4 

•it r 
• • i 


Vik k 

Xt! r 

lUtk . . -* 

• ii i i 

v*j j 

rrJ-JJ- • ■ "» 
i •• i • 

• Ml f 

• in • 




1 1 Dill * 1 I I • 1 III 

» * • t I 1 ■ 1 It 

1 • » • Ml • • ••ttftl 


t • • 

i « • 






mi i 



FIG 15 
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P/N diagram of the target filter 
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